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Abstract — The  Wigner-Ville  Hough  transform  (WVHT)  has 
been  applied  to  detect  and  estimate  the  parameters  of  linear 
frequency-modulated  continuous-wave  (LFMCW)  low  probabil¬ 
ity  of  intercept  (LPI)  radar  waveforms.  The  WVHT,  which  is 
optimal  for  a  single  linear  frequency  modulated  (LFM)  signal, 
becomes  sub-optimal  when  applied  to  LFMCW  signals  since  the 
observed  waveform  is  composed  of  concatenated  LFM  pulses. 
We  formulate  the  detection  and  estimation  problem  to  take  into 
account  the  multiple  pulses  that  are  available  in  an  observation 
interval  at  the  intercept  receiver.  The  new  algorithm,  called  the 
periodic  WVHT  (PWVHT),  is  shown  to  significantly  outperform 
the  WVHT  for  LFMCW  signals. 

I.  Introduction 

In  electronic  warfare,  knowledge  of  enemies’  electronic 
capabilities  is  desired.  Electronic  warfare  includes  the  inter¬ 
ception  and  analysis  of  radar  signals  (ELINT)  to  determine 
information  about  their  capabilities  [1].  The  general  process  of 
ELINT  includes  a  detection  and  estimation  stage  in  which  the 
sensor  decides  whether  signals  of  interests  (SOIs)  are  present 
and,  if  they  are,  estimates  their  parameters,  such  as  the  carrier 
frequencies  and  pulse  durations.  The  main  topic  of  this  paper 
is  the  detection  and  parameter  estimation  of  linear-frequency- 
modulated  continuous  wave  (LEMCW)  radar  signals. 

Various  signal-processing  techniques  have  been  investigated 
for  the  detection  and  estimation  of  low  probability  of  intercept 
(LPI)  radar  signals  [2]-[6],  including  LEMCW  signals.  (Eor  a 
detailed  discussion  of  LPI  radar  signals,  see  [7]).  In  general, 
the  detection  and  estimation  of  linear  PM  (LEM)  signals  has 
been  well  studied  in  the  literature.  One  of  the  most  prominent 
techniques  is  the  Wigner-Ville  Hough  Transform  [8],  [9].  The 
Wigner-Ville  Hough  transform  (WVHT),  the  Hough  transform 
(HT)  of  the  Wigner-Ville  distribution  (WVD),  has  been  shown 
to  be  equivalent  to  the  generalized  likelihood  ratio  test  (GLRT) 
and  maximum  likelihood  estimator  (MLE)  (i.e.,  it  is  asymptot¬ 
ically  optimal)  in  detection  and  estimation  of  an  LEM  signal 
[10].  A  sawtooth  LEMCW  signal  is  a  periodic  extension  of  an 
LEM  signals,  i.e.,  the  frequency  within  a  modulation  period 
of  an  LEMCW  radar  signal  changes  linearly,  as  illustrated  in 
Pigure  1.  The  WVHT  is  not  optimal  for  the  detection  and 
estimation  of  LEMCW  radar  signals  since  they  are  not  LEM 
when  we  consider  the  whole  observation  time. 

In  radar  signal  interception,  the  observation  time  can  be 
controlled  to  the  advantage  of  the  intercept  receiver:  by  having 
a  longer  observation  time,  more  pulses  are  available  within 


Fig.  1.  Available  LFMCW  radar  signal  within  an  observation  time  Tohs 

the  observation  signal  and  thus  there  is  more  energy  available 
for  detection.  Intuitively,  the  detection  and  estimation  perfor¬ 
mance  ought  to  improve  as  the  observation  time  is  increased. 
Unfortunately,  this  is  not  the  case  with  the  WVHT,  since  it  is 
tailored  for  LEM  signals. 

In  this  paper  we  present  a  detection  and  estimation  tech¬ 
nique  named  the  periodic  WVHT  (PWVHT)  that  is  optimal  for 
LEMCW  signals.  The  algorithm  is  similar  to  the  WVHT  ex¬ 
cept  that  it  searches  for  patterns,  such  as  a  sawtooth  waveform, 
in  the  time-frequency  (TP)  image.  It  is  thus  able  to  accumulate 
signal  energy  over  the  entire  observation  interval,  unlike  the 
traditional  WVHT.  The  PWVHT  is  shown  to  outperform  the 
WVHT  applied  to  a  LEMCW  signal. 

H.  Interceptor  eor  LEMCW  Signals 

We  assume  that  the  signal  processing  block  in  the  intercep¬ 
tor  is  presented  with  a  noisy  discrete-time  observation  signal 
r[n],  n  =  0, 1, ...,  V  —  1,  where  N  is  the  number  of  samples. 
If  SOIs  are  present,  these  bandpass  continuous  LEMCW  radar 
signals  are  denoted  as  s' {t).  Let  s{t)  be  the  signal  s'{t)  down- 
converted  to  a  known  intermediate  frequency,  which  is  then 
sampled  with  sampling  period  A  (sec). 

The  LEMCW  radar  signal  can  be  written  as 

s{t)  =  ^ei(¥^+27r/it+7r^mod(t+r,T)2)^ 

where  A  is  the  amplitude,  fi  (Hz)  is  the  initial  frequency  after 
down-conversion,  and  g  =  ^  (^)  is  the  chirp-rate,  where  B 
(Hz)  is  the  bandwidth,  r  (sec)  is  an  initial  time-offset,  and 
T  (sec)  is  the  modulation  period.  The  symbol  ip  represents 
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the  initial  phase,  which  is  distributed  uniformly  f/(— 7r,7r), 
and  mod(t,  T)  is  a  shorthand  notation  for  (t  modulo  T).  The 
observation  time  t  (sec)  at  the  intercept  receiver  is  bounded  by 
Tobs,  i-e.,  0  <  t  <  Tobs-  A  time  vs.  frequency  (t,  /)  illustration 
of  the  LFMCW  radar  signal  s{t),  along  with  its  parameters  is 
given  in  Figure  1. 

The  discrete-time  LFMCW  radar  signals  that  may  be  present 
in  r[n]  can  be  written  as  5[n]  =  s{t)  \t=nA-  Given  a  noisy 
observation  signal  r[n],  the  decision  problem  is  to  decide 
whether  hypothesis  Hq  or  Hi  is  true,  where 

Ho  :r[n]  =  w[n]  (2) 

Hi  :  r[n]  =  5[n]  w[n]  (3) 

The  noise  w  [n]  is  zero-mean,  stationary,  complex  white  Gaus¬ 
sian,  i.e.,  A/’(0,cr^).  If  the  decision  is  that  hypothesis  Hi 
is  true,  then  an  estimate  of  the  true  parameter  vector  Q  = 
(fi^g^T^r)  is  desired.  The  amplitudes  and  phases  of  the 
signals  are  considered  nuisance  parameters. 

III.  Periodic  Wigner-Ville  Hough  Transform 

The  crucial  drawback  of  the  WVHT  when  used  for  the  de¬ 
tection  and  estimation  of  LFMCW  radar  signals  at  the  intercept 
receiver  is  that  there  no  longer  is  gain  in  performance  when 
the  observation  time  is  increased.  In  this  section  we  introduce 
the  periodic  Wigner-Ville  Hough  transform  (PWVHT),  which 
is  based  on  the  WVHT  but  tailored  to  LFMCW  signals.  The 
transformation  from  the  time  and  frequency  (t,  /)  domain 
to  the  initial  frequency  and  chirp  rate  (/,^)  domain  in  the 
WVHT  is  done  by  summing  over  straight  lines,  parameterized 
by  (/,^),  in  the  WVD.  For  an  arbitrary  signal  x[n]  it  can  be 
written  as 

oo 

W,[n,f)  =  2  Y,  C,,[n,k]F{f,k],  (4) 

k=  —  oo 

where  F(/,  k]  =  g-i^Tr/ZcA  matching  function, 

Gccy  [n,  k]  =  [n—k]  is  the  instantaneous  discrete  time 

cross-correlation,  and  the  symbol  stands  for  the  complex 
conjugate.  The  WVD  is  basically  the  Fourier  transform  of 
the  instantaneous  auto-correlation  function  with  respect  to  k. 
For  a  single  LFM  signal  5[n]  =  ^  \t=nA, 

the  instantaneous  auto-correlation  function  is  Css[n,k]  = 
^2gi47r/cA(/,+^nA)  the  WVD  of  an  LFM  signal  is 

oo 

Ws[n,f)=  Y  ^2e-i47rfeA(/-/i-gnA)^  (5^ 

k=  —  oo 

Notice  that  each  term  in  the  WVD  is  maximized  only  when 
it  is  evaluated  at  the  signal  parameters  f  =  fi  gnA  at  any 
given  time  n  since  the  exponential  terms  cancel  out. 

The  WVHT  can  be  written  as 

N/2-1  n 

WH, if,  9]=  Y 

n=0  k=—n 

N—1  N—l—n 

+  E  E  Cxx[n,k]F{f  +  gnA,k].  (6) 

n=N/2  k=-(N-l-n) 


Therefore,  the  discrete  WVHT  of  a  discrete  LFM  signal  is 

WH,  [f,  5]  =  E  E  (7) 

n  k 

and  is  maximized  only  when  evaluated  at  f  =  fi  and  g  =  g, 
which  results  in 

WH,[fi,g]  =  ^-.  (8) 

It  can  be  seen  that  the  canceling  of  complex  exponentials  in 
(5)  and  (7)  when  evaluated  at  the  signal  parameters  is  a  key 
in  detection  and  estimation  of  an  LFM  signal  by  the  WVHT. 

In  the  case  of  an  LFM  signal,  F{fA]  and  Cxx[f^ik]  are 
a  good  match.  The  Hough  transformation  involves  integrating 
the  WVD  along  the  frequency  law  pattern  of  the  signal  in  the 
TF  plane  for  all  the  possible  values  for  the  parameters,  and 
only  at  the  correct  signal  parameters  is  the  sum  maximized, 
which  is  equivalent  to  when  F{f,k]  completely  cancels  out 
the  exponential  term  in  Cxx[^^  k]. 

The  WVHT  is  capable  of  detecting  and  estimating  LFM 
signals  even  when  more  than  one  signal  is  present,  as  seen 
in  Figure  2  where  three  LFM  signals  are  processed,  yet  it 
is  not  optimal  for  more  than  one  LFM  signal.  The  main 
reason  is  that  the  matching  function  cannot  completely  remove 
the  exponential  term  in  Cxx[riA]-  For  instance,  the  discrete 
instantaneous  auto-correlation  of  a  noise-free  observation  r[n] 
consisting  of  two  LFM  signals  si[n]  and  52 [n]  is 

Crr[n,  k]  =  Cs,s,  [n,  k]  +  [n,  k] 

+CsiS2  [n,  k]  +  Cs2Si  [n,  k]  (9) 

and  it  is  easy  to  see  that  F{f,k]\j:^j:_^_-^  cannot  completely 
eliminate  the  exponential  terms  in  (9).  When  the  value  of  / 
in  F(/,  k]  matches  one  of  the  signal  parameters,  it  cancels 
out  one  of  the  exponential  term,  either  in  Cs^sii^A]  or 
Cs2S2  k],  yielding  prominent  peaks  in  (6).  It  also  shows  why 
it  is  sub-optimal  since  F{f,k]  can  only  cancel  one  exponential 
term  at  a  time.  Thus,  the  WVHT  of  M  LFM  signals  yields  M 
peaks,  which  need  to  be  searched  by  a  peak  finding  algorithm. 

The  WVHT  applied  to  LFMCW  signals  is  also  sub-optimal 
because  the  matching  function  F{fA]  cannot  completely 
cancel  out  the  exponential  terms  in  (6)  when  multiplied  with 
the  instantaneous  auto-correlation  function  Cxx[i^A]  of  an 
LFMCW  radar  signal,  i.e.,  they  are  not  a  good  match.  For 
an  LFMCW  radar  signal 

s[n]  =  ^^jA-t^^fi(nA)-\-7Tgmod(nA-\-T,T)^)^ 

The  instantaneous  auto-correlation  function  is 

Css[n,  k]  =  ex.p{j7r{4fikA 
-h^[mod((n  +  /c)A+r,  T)^  —  mod((n  — /c)A+r,  T)^])}(11) 

Thus,  the  WVHT  in  (6)  results  in 

WHs[f,g]  =  ^y]^2exp{j7r(4(/j-/-5n)A:A 

n  k 

+^[mod((n  +  k)A  +  r,  T)^  —  mod((n  —  k)A  +  r,  T)^])}, 
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Fig.  2.  WVD  and  WVHT  parameterized  by  (/,  g)  of  a  single  LFMCW  signal  formed  by  three  LFM  pulses.  Noise-free:  (a)  WVD,  (b)  WVHT.  Noisy  with 
SNR^  =  —10  dB:  (c)  WVD,  (d)  WVHT.  In  this  and  subsequent  figures,  the  frequency  is  normalized  by  1/A,  the  modulation  period  and  time  offset  by 
and  the  chirp-rate  by  the  maximum  chirp-rate  searched. 


and  the  exponentials  do  not  cancel  out  for  any  values  of  /  and 
g.  Thus,  the  summation  results  in  a  smaller  peak. 


The  optimal  algorithm  for  one  LFMCW  signal  uses  a  match¬ 
ing  function  that  can  completely  cancel  out  the  exponential 
terms  when  multiplied  with  (11).  We  define  the  PWVHT  of 
an  arbitrary  complex  sequence  x[n],  n  =  0, 1, ...,  —  1  to  be 

N/2-1  n 

z.pi  =  E  s;  k]F n,,/c[^] 

n=0  k=—n 
N-1  N-l-n 

E  E  Cxx[n,k]Fn,kM,  (12) 

n=N/2  k=-{N-l-n) 

where  ft  =  (/,^,T,f).  We  assume  that  N  is  even  and  the 


matching  function  is 

Fn,k[^]  =  exp{-ji'7r(4/fcA  +  5[inod((n  +  k)A  +  - 

mod((n  - /i:)A  +  f,  f)^])},  (13) 

which,  when  evaluated  at  the  true  signal  parameters  Q  =  ft  = 
(/i,  T,  r),  is  the  complex  conjugate  of  the  exponential  term 
of  (1 1).  The  PWVHT  is  optimal  as  it  is  equivalent  to  the  GLRT 
and  the  estimates  it  yields  are  MLE. 

For  an  LFMCW  radar  signal  at  the  intercept  receiver,  the 
peak  of  the  PWVHT  is  achieved  at  the  signal  parameters  ft. 
For  instance,  if  we  substitute  the  LFMCW  signal  in  (10)  into 
(12),  and  evaluating  it  at  the  true  signal  parameters,  it  yields 
Zs[ft]  =  ^ ^  ,  which  is  the  maximum  and  equal  to  (8).  For 
other  values,  i.e.,  ft  ^  ft,  the  matching  function  does  not 
cancel  out  the  exponential  terms,  resulting  in  a  smaller  value. 
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Fig.  3.  Slices  of  Zr[^]  of  a  noisy  LFMCW  signal  in  Figure  2:  (a)  slice  at 
the  true  modulation  period  and  time  offset  of  the  signal,  (b)  slice  at  the  true 
initial  frequency  and  chirp  rate  of  the  signal.  SNR^  =  —10  dB. 


IV.  Simulation  Results  and  Analysis 


To  illustrate  the  performance  of  the  PWVHT,  Z[^],  we 
apply  it  to  the  simulated  LFMCW  radar  signal  shown  in 
Figure  2.  When  the  PWVHT  of  the  noisy  signal  is  evaluated, 
a  prominent  peak  appears,  as  seen  in  Figure  3,  where  slices 
of  the  Z[^]  are  shown.  The  right  half  of  Figure  3(b)  is 
smooth  because  only  time  offsets  less  than  the  modulation 
periods  need  to  be  searched.  Each  plot  indicates  a  peak, 
revealing  the  correct  number  of  signals  (one  instead  of  three). 
In  addition,  the  location  of  the  peak  gives  an  estimate  of  the 
signal  parameters.  The  detection  and  parameter  estimation  of 
the  LFMCW  signals  can  then  be  achieved  by  searching  for 
peaks  in  the  (l  domain. 
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Fig.  4.  ROC  plot  of  the  PWVHT  and  the  WVHT  for  a  single  LFMCW 
signal,  using  N  =  330  and  Tm  =  IlOA  (3  concatenated  LFW  signals) 


The  output  SNR  can  be  canculated  to  be 

.  2 


SNRpwVHT — 


(«■ 


N^  +  1 


YSNR^ 
NSNRi  + 1  ■ 


(14) 


This  is  a  remarkable  result  because  when  the  input  SNR  is 
high,  (SNRi  =  ^  >>  1),  the  output  SNR  is  approximately 

SNRpwvHT  =  ,  which  indicates  that  it  is  almost  as 

good  as  a  matched  filter’s  output  SNR  (SNRmf  =  A^SNR^), 
but  with  approximately  3  dB  loss  due  to  quadratic  detection. 
Most  importantly,  when  the  input  SNR  is  very  low  (SNR^  << 
1),  detection  can  still  be  accomplished  by  increasing  number 
of  samples  N  and  benefiting  from  processing  gain.  Not 
surprisingly,  the  PWVHT  is  equivalent  to  the  GLRT  and  the 
MLE  for  an  LFMCW  radar  signal  at  the  intercept  receiver. 
The  detection  performance  of  the  PWVHT  is  shown  using  a 
received  operating  characteristic  (ROC)  plot  in  Figure  4.  The 
probability  of  detecting  a  simulated  LFMCW  signal  using  the 
PWVHT  is  compared  to  the  same  metric  using  a  standard 
WVHT  for  two  input  SNR  values,  as  the  probability  of  false 
alarm  varies.  The  LFMCW  algorithm  is  clearly  superior  to 
the  WVHT,  as  predicted  since  the  signal  is  LFMCW  and  not 
purely  LFM. 


Although  the  PWVHT  is  superior  to  the  WVHT  for  detec¬ 
tion  and  estimation  of  a  signal  that  is  composed  of  multiple 
LFM  signals,  the  improved  performance  comes  at  a  cost 
in  computational  complexity.  However,  for  LPI  application, 
the  PWVHT  is  realizable.  For  detection  and  estimation  of 
LFMCW  radar  signal,  the  search  is  usually  bounded,  i.e., 
there  is  a  priori  knowledge  about  the  range  of  possible  initial 
frequencies,  chirp-rates,  and  modulation  periods,  which  allows 
one  to  restrict  the  search  in  the  parameter  space  and  reduce 
the  amount  of  computations. 
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V.  Conclusion 

In  this  paper  we  present  a  new  signal  processing  technique, 
the  PWVHT,  for  detection  and  estimation  of  a  LFMCW  radar 
signal.  Within  an  observation  session  at  an  intercept  receiver, 
the  received  signal  is  composed  of  multiple  LFM  pulses  if  an 
SOI  is  present.  Currently  available  techniques  are  suboptimal 
since  they  are  tailored  for  the  detection  and  estimation  of  a 
single  pulse.  The  PWVHT  technique  is  equivalent  to  the  GLRT 
and  MLE  for  an  LFMCW  radar  signal,  and  thus  provides  an 
optimal  solution.  The  improvement  in  performance  comes  at 
a  cost  in  computational  complexity. 

The  idea  in  the  PWVHT  can  be  extended  to  other  signal 
types  (i.e.,  other  signals  composed  of  multiple  LFM  pulses) 
by  designing  a  matching  function  that  can  completely  remove 
the  exponential  term  when  multiplied  with  its  instantaneous 
auto-correlation  function.  This  warrants  further  research. 
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